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Abstract. 

We discuss three aspects of the interaction between GRB and their surroundings. The 
illumination of the progenitor remnant and/or the surroundings by the X-ray afterglow 
continuum can produce substantial Fe K-alpha line and edge emission, with implica- 
tions for the progenitor model. The presence of large dust column densities, capable 
of obscuring the GRB optical afterglow, will lead to characteristic delayed X-ray and 
far-IR light curve signatures. Pair production induced by the initial gamma-rays in the 
nearby environment will modify the initial spectrum and the afterglow light curve, and 
the magnitude of these changes provides a diagnostic for the external density. 
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I FE X-RAY LINES FROM GRB PROGENITORS 

Important clues for identifying the nature of the progenitors of the long (t>2 s) GRBs may be 
available from the recent report at a 4.7c level of X-ray Fe line features in the afterglow after 1.5 
days of the gamma-ray burst GRB 991216 [10], as well as similar detections at the 3a level in 5 
other bursts with Beppo-SAX and ASCA. X-ray atomic edges and resonance absorption lines are 
theoretically expected to be detectable from the gas in the immediate environment of the GRB, 
and in particular from the remnants of a massive progenitor stellar system [15,14,16]. 

A straightforward interpretation [10] of the GRB 991216 observation would imply a mass >,0.1 — 
1M Q of Fe at a distance of about 1-2 light-days, possibly due to a remnant of an explosive event 
or supernova which occurred days or weeks prior to the gamma-ray burst itself (a 'supranova', 
[10,12]). The long time delay is necessary both to get the relatively massive, slow moving ejecta 
out to few light-day distances (to explain the line appearance at a few days with light travel 
arguments), and in order to get the initial Ni and Co to decay to Fe (~ 55 days). This requires 
a two-step process, in which an initial supernova leads to a temporarily stabilized neutron star 
remnant, which after weeks collapses to a black hole leading to a canonical burst ( [11,12]). It is 
unclear whether fall-back from the supernova leading to the second collapse to a BH could occur 
with such a (~ weeks) long delay (e.g. [18]). Another possibility is that a massive progenitor has 
previously emitted a copious wind (A/ i >lO~ 4 Af /yr), which would need to be unusally Fe-rich 
and highly inhomogeneous ( [14]; c.f. [10]). 

An alternative, and perhaps less restrictive scenario for such Fe lines [17] involves an extended, 
possibly magnetically dominated wind from a GRB impacting the expanding envelope of a massive 
progenitor star. This could be due either to a spinning-down millisecond super-pulsar or to a 



highly-magnetised torus around a black hole (e.g. [13]), which could produce a luminosity that 
was still, one day after the original explosion, as high as L m ~ W 47 t^ y 3 ergs. An outflow with 
such a dependence can also be powered by accretion of fall-back material onto a central black hole 
[18]. This jet luminosity may not dominate the continuum afterglow; but its impact on the outer 
portions of the expanding stellar envelope at distances ,<10 13 cm, even with just solar abundances, 
can be efficiently reprocessed into an Fe line luminosity comparable to the observed value, together 
with a contribution to the X-ray continuum. Under this interpretation, the dominant continuum 
flux in the afterglow, even in the X-ray band, is still attributable to a standard decelerating blast 
wave. 

The relativistic magnetised wind from the compact remnant would develop a stand-off shock 
before encountering the envelope material, and shocked relativistic plasma would be deflected 
along the funnel walls. Non-thermal electrons will be accelerated behind the standoff shock in the 
jet material; the transverse magnetic field strength (which decreases as 1/r in an outflowing wind) 
would be of order 10 4 G at 10 13 cm - strong enough to ensure that the shock-accelerated electrons 
cool promptly, yielding a power- law continuum extending into the X-ray band. Some of these 
X-rays would escape along the funnel, but at least half (the exact proportion depending on the 
geometry and flow pattern) would irradiate the material in the stellar envelope. Pressure balance 
in the shoked envelope wall implies densities of n e = aL m /6TTr 2 ckT ~ 10 17 a L 47 r^T^ 1 cm~ 3 , 
where a ~ 1 is a geometric factor, and the recombination time for hydrogenic Fe in the funnel walls 
photoionized by the non-thermal continuum is t rec = 6 x lO~ 6 T^ 2 n^ ~ 6 x 10^ 6 aL~ n L i7 rl 3 T^ 2 s. 
Standard calculations of photoionization of optically-thin slabs (e.g. [7]) show that the equivalent 
width of the Fe K-alpha line, for solar abundances, is about 0.5 kev, or twice as strong if the Fe 
has ten times solar abundances. These results are applicable provided that the ionizing photons 
encounter a Fe ion before being scattered by free electrons i.e. provided that tt = GT&%n e <\. 
Under these conditions the Fe K-a photon flux is about 0.1 of the X-ray continuum [17], N~LFe ~ 
10 54 L47/3 ph/s, where f3 < 1 is the ratio of ionizing to MHD luminosity. This line luminosity 
compares well with Fe line luminosity 6 x 10 52 ph/s observed t ~ 1.5 day after the GRB 991216 
burst by [10]. 

The total amount of Fe needed to explain the observed K-a line flux, arising in a thin layer of 
the funnel walls of a collapsar model, amounts to a very modest mass of Mp e ~ 10~ 8 Mq, which 
could be Fe synthesized in the core. The Fe-enriched core material can easily reach a distance 
comparable to r ~ 10 13 cm in 1 day for an expansion velocity below the limit v ~ 10 9 cm s _1 
inferred by [10] from the line widths. Even without this, a solar abundance (10 _5 Mq of Fe) in the 
envelope is sufficient to explain the observations. The initial, energetic portion of the relativistic 
jet, with a typical burst duration of 1 — 10 s, will rapidly expand beyond the stellar envelope, 
leading in the usual way to shocks and a decelerating blast wave. A continually decreasing fraction 
of energy, such as put out by a decaying magnetar, may continue being emitted for periods of a 
day or longer, and its reprocessing by the stellar envelope can be responsible for the observed Fe 
line emission in GRB 991216. Since the energy in this tail can decay faster than t -1 , the usual 
standard shock gamma-ray and afterglow scenario need not be affected, being determined by the 
first 1-10 s worth of the energy input. 

II DUSTY GRB DELAYED XR/IR AFTERGLOWS 

For GRB in large star forming regions, a significant fraction of the prompt X-ray emission 
will be scattered by dust grains. Since dust grains scatter X-rays by a small angle, time delays 
of the scattered x-rays will be small (minutes to days, depending on the X-ray energy and the 
grain size). If the dust column density is substantial, the softer part of the X-ray afterglow on 
the above timescales will be dominated by the dust scattering, the direct X-ray emission from the 



blast wave being weaker. This intermediate time, soft(cr) X-ray light curve will be steeper than 
the unscattered X-ray afterglow. 

As a specific example [6], consider a typical GRB whose unscattered X-ray light curve is 
parametrized as F (t) = [l + (t/100s)]/[l + (t/100s) 2,3 ], with an arbitrary normalization depending 
on the X-ray energy band. This is represented by the thin line in Fig.l. We assume that the 
GRB occurs in a large star forming region, of typical radius R about lOOpc, where the dust grain 
populations and optical depths are close to what is observed in our Galactic center region. Thus 
for numerical estimates we assume that (1) visual extinction is ~ 10, (2) X-rays are scattered 
preferentially by those dust grains whose size is in the range a ~ 0.06/im, (3) the optical depth 

to dust scattering at the X-ray energy e is r(e) = 3 ( lk ^ v ) . At X-ray optical depths less than 
few, dust grains of size a will scatter X-rays of energy e by an angle 9 ~ 0.2A/a, where A is 

the X-ray wavelength, 6(e) - 4x 1CT 3 (q^Ts) (iiSv) _1 - The time la § is 1 ~ R02 / 2c ' or 

*oo - 9 >< iq4s (ofc) 2 (i^r 2 (to) • 

At 2 keV, the optical depth is r ~ 1. The time lag is t ~ 2 x 10 4 s. The scattered flux is 
F s ~ rf/t ~ 0.03. The unscattered flux at 2 x 10 4 s is F ~ 10~ 3 . In the time interval from hours 
to weeks, the dust scattering dominates the afterglow, and, as shown in Fig.l, the afterglow is 
approximately a power law F oc t~ 175 [6]. This is because dust grains of radius a < 0.06/im will 
scatter the prompt emission with longer time lags, t oc a~ 2 , and with smaller optical depths r. 
To calculate r, we take a standard dust grain size distribution where the number of grains of size 
of order a is oc a~ 2 - 5 For a scattering cross section oc a 4 the optical depth is r oc a 15 oc t -0 - 75 , so 
the flux F oc t' 1 - 75 . 




FIGURE 1. Dust-scattered X-ray af- 
terglow. Thin line: unscattered X-ray 
flux. Thin dashed line: scattered X-ray 
flux. Thick line: total flux. The flux 
normalization is arbitrary, the relative 
fluxes correspond to the example dis- 
cussed in the text for an energy of 2 
keV [6]. 



A GRB in such a highly obscured star-forming region should lead to specific signatures in the 
X-ray afterglow, i.e. a bump in the X-ray light curve at energies e ~ 2 — 3 keV, hours to days 
after the burst [6] . This X-ray signature is expected for bursts which do not produce a detectable 
optical transient (OT). Such OT-less, X-ray peculiar GRBs will also lead to thermal reemission 
and scattering of the O /UV flux causing a delayed IR emission, as is the case also for partially 
absorbed bursts [8,9]. For an isotropic total burst energy E ~ 10 53 erg at a redshift z ~ 1 the 
normalization of the X-ray flux for the burst of Figure 1 would be F x ~ 10 -9 erg cm -2 s _1 
keV -1 for t < 100 s, in the usual range of X-ray afterglow fluxes detected by Beppo-SAX. The dust 
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reradiation occurs beyond the sublimation radius R s ~ 10 L 4 g pc at wavelengths A>2(1 + z)/im, 
where 10 49 L4g erg/s is the early UV component of burst afterglow [8]. The time delay associated 



with the reradiated flux is tiR ~ (R s /2c)6j where 6j = 10 _1 ^_i. At z ~ 1 the corresponding 

infrared flux at 2.2 fim would be F 2 .2^m ~ L^Oj /[AnDi i (R a /2c)9ji>] ~ O.SL^ 2 /uJy, independent 

of 6*j, or mjf ~ 23.3 compared to Vega [6], approximately constant for a time i/^ ~ 5 x 10 6 6 2 _ 1 L\g 2 
s. Such 7-ray detected GRBs with anomalous X-ray afterglow behavior and no OT may be used 
as tracers of massive stellar collapses. It may thus be possible to detect star-forming regions out 
to rcdshifts larger than achievable with O/IR techniques, since typical GRB 7-ray, X-ray and IR 
fluxes can in principle be measured out to z ~ 10 — 15. 



Ill PAIR PRODUCTION IN GRB ENVIRONMENTS 

Gamma-ray burst sources with a high luminosity can produce pair cascades in their envi- 
ronment as a result of back-scattering of a seed fraction of their original hard spectrum. New 
pairs can be made as some of the initial energetic photons are backscattered and interact with 
other incoming photons. Previous work on this investigated the acceleration of new pairs for 
a particular fireball model [1,4], the effect of pair formation for a low compactness parameter 
external shock model of GRB [3], and Compton echos produced by pairs [2]. Here we discuss a 
simplified analytical treatment [5] of pair effects from 7-rays arising in internal shocks in a wind; 
the remaining wind energy drives a blast wave which decelerates as it sweeps up the external 
medium, and gives rise to the afterglow emission. The 7-rays would propagate ahead of the blast 
wave, leading to pair production (and an associated deposition of momentum) into the external 
medium. The pair cascades saturate after the external (pair-enriched) medium reaches a critical 
bulk Lorentz factor, which is generally below that of the original relativistic wind. For external 
baryonic densities similar to those in molecular clouds the pairs can achieve scattering optical 
depths t±<1. Even for less extreme external densities the effect of the additional pairs can be 
substantial, increasing the radiative efficiency of the blast wave and leading to distortions of the 
original spectrum. This provides a potential tool for diagnosing the compactness parameter of 
the bursts and thus the radial distance at which shocks can occur. It also provides a tool for 
diagnosing the baryonic density of the external environment, and testing the association with 
star- forming regions. 

For the maximum Lorentz factor to which an can be accelerated by scattering, and the 
maximum Lorentz factor at which back-scattered photons can still make new pairs, one finds 
two regimes defined by the effective duration of the light pulse seen by the screen of accel- 
erated pairs. At low radii (wind regime) the effective duration is the burst duration t w ; for 
large radii (impulsive regime), the effective duration is At ~ rjcT\. For an incident photon 
number index (3 = 2, in the former T± oc r~ 1//3 and in the latter T± cx r~ 2 . The critical 
radius and Lorentz factor for the transition between the wind and the impulsive regimes are 
[5] r c = 5 x 10 14 L^ // 5 5 ^ / ' 1 5 , F c = 3 x lO 1 !,^^"^ 5 . The maximum radius at which pair 
cascades cut off is re ~ (Ar*ct w /i) 1 / 2 ~ 4 x ^ 15 ^]ih6^wi cm. Before the pairs start acceler- 
ating, assuming they are held back by the environmental protons through magnetic fields, an 
initial cascade amplification fator k p ~ (m p /m e ) is achieved. After the mean mass per scat- 
terer drops to a value comparable to the electron mass, before reaching re a further amplification 
k a ~ 2 s ~ 50 (where s ~ logT c / \og2) is possible, so the total pair amplification factor is [5] 
k c = k p k a (r c ) ~ (m p /m e )T c ~5x lO 4 /^^^ 5 . The maximum pair optical depth at r c , which 
is prevented from exceeding t± ~ 1 by self-shielding, is achieved for external densities n p >n p , c , 

where n p , c ~ 10 5 L~5^ 5 ir,i /5 cnr 3 . 

The external density and the initial Lorentz factor r\ determine when the outer shock and the 
reverse shock become important and whether this happens within the radius already polluted 
with pairs. If i]<ri/ct w the external shock responsible for the afterglow occurs beyond the region 



"polluted" by new pairs, and otherwise the afterglow shock may experience, after starting out 
in the canonical manner, a "resurgence" or second kick as its radiative efficiency is boosted by 
running into an e ± -enriched gas [5]. 

Additional effects are expected when r± — ► 1, for external baryon density n p >n c ^ p ~ 

^Q 5 L~5Q 5 t~^ 5 cm~ 3 at radii r < re. Such high densities could be expected if the burst is as- 
sociated with a massive star in which prior mass loss led to a dense circumstellar envelope. The 
pair optical depth saturates to r± ~ 1 and in addition to an increased efficiency and softer spec- 
trum of the afterglow reverse shock, the original gamma-ray spectrum of the GRB will be modified 
as well. One of the consequences of such a critical external density leading to t± ~ 1 would be 
the presence of an X-ray quasi-thermal pulse, whose total energy may be a few percent of the 
total burst energy [5]. 
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